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We consider f(R,T) modified theories of gravity, where the gravitational Lagrangian is given 
by an arbitrary function of the Ricci scalar R and of the trace of the stress-energy tensor T. We 
obtain the gravitational field equations in the metric formalism, as well as the equations of motion 
for test particles, which follow from the covariant divergence of the stress-energy tensor. Generally, 
the gravitational field equations depend on the nature of the matter source. The field equations 
of several particular models, corresponding to some explicit forms of the function f(R,T), are also 
presented. An important case, which is analyzed in detail, is represented by scalar field models. We 
write down the action and briefly consider the cosmological implications of the / (J?, T"^) models, 
where T^ is the trace of the stress-energy tensor of a self-interacting scalar field. The equations of 
motion of the test particles are also obtained from a variational principle. The motion of massive 
test particles is non-geodesic, and takes place in the presence of an extra force orthogonal to the 
four-velocity. The Newtonian limit of the equation of motion is further analyzed. Finally, we provide 
a constraint on the magnitude of the extra-acceleration by analyzing the perihelion precession of 
the planet Mercury in the framework of the present model. 

PACS numbers: 04.20.Cv, 04.50.Kd, 98.80.Jk, 98.80.Bp 



I. INTRODUCTION 

The recent observational data (T], on the late-time acceleration of the Universe and the existence of dark matter 
have posed a fundamental theoretical challenge to gravitational theories. One possibility in explaining the observations 
is by assuming that at large scales the Einstein gravity model of general relativity breaks down, and a more general 
action describes the gravitational field. Theoretical models, in which the standard Einstein-Hilbert action is replaced 
by an arbitrary function of the Ricci scalar R jsj , have been extensively investigated lately. The presence of a late-time 
cosmic acceleration of the Universe can indeed be explained by f(R) gravity The conditions of the existence of 
viable cosmological models have been found in [j| , and severe weak field constraints obtained from the classical tests 
of general relativity for the Solar System regime seem to rule out most of the models proposed so far [y, [7[. However, 
viable models, passing Solar System tests, can be constructed pH-tlH- f{R) models that satisfy local tests and unify 
inflation with dark energy were considered in 12]. In the framework of f(R) gravity models the possibility that the 
gala ctic dynamic of massive test particles can be understood without the need for dark matter was considered in 
|13l - [T7| . For reviews of f(R) generalized gravity models see [!, Ojj]. 
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A generalization of f(R) modified theories of gravity was proposed in fl9| . by including in the theory an explicit 
coupling of an arbitrary function of the Ricci scalar R with the matter Lagrangian density L m . As a result of the 
coupling the motion of the massive particles is non-geodesic, and an extra force, orthogonal to the four- velocity, arises. 
The connections with Modified Newtonian Dynamics (MOND) and the Pioneer anomaly were also explored. This 
model was extended to the case of the arbitrary couplings in both geometry and matter in [20| . The astrophysical 
and cosmological implications of the non-minimal coupling matter-geometry coupling were extensively investigated in 
[2ll |22| , and the Palatini formulation of the non- minimal geometry-coupling models was considered in [23j j . In this 
context, a maximal extension of the Hilbert-Einstein action was proposed in [24|, by assuming that the gravitational 
Lagrangian is given by an arbitrary function of the Ricci scalar R and of the matter Lagrangian L m . The gravitational 
field equations have been obtained in the metric formalism, as well as the equations of motion for test particles, which 
follow from the covariant divergence of the stress-energy tensor. 

A specific application of the latter f(R, L m ) gravity was proposed in (25j . which may be considered a relativistically 
covariant model of interacting dark energy, based on the principle of least action. The cosmological constant in 
the gravitational Lagrangian is a function of the trace of the stress-energy tensor, and consequently the model was 
denoted "A(T) gravity" . It was argued that recent cosmological data favor a variable cosmological constant, which 
are consistent with A(T) gravity, without the need to specify an exact form of the function A(T) [25[. A(T) gravity 
is more general than the Palatini f{R) gravity, and reduces to the latter when we neglect the pressure of the matter. 

It is the purpose of the present paper to consider another extension of standard general relativity, the f(R,T) 
modified theories of gravity, where the gravitational Lagrangian is given by an arbitrary function of the Ricci scalar R 
and of the trace of the stress-energy tensor T. Note that the dependence from T may be induced by exotic imperfect 
fluids or quantum effects (conformal anomaly). As a first step in our study we derive the field equations of the 
model from a variational, Hilbert-Einstein type, principle. The covariant divergence of the stress-energy tensor is also 
obtained. The f(R, T) gravity model depends on a source term, representing the variation of the matter stress-energy 
tensor with respect to the metric. A general expression for this source term is obtained as a function of the matter 
Lagrangian L m . Therefore each choice of L m would generate a specific set of field equations. Some particular models, 
corresponding to specific choices of the function f(R,T) are also presented, and their properties are briefly discussed. 
In fact, we also demonstrate the possibility of reconstruction of arbitrary FRW cosmologies by an appropriate choice 
of a function f(T). Scalar fields play a fundamental role in cosmology, i.e., as possible explanations for inflation, late 
time acceleration, or dark matter, respectively. Therefore, we introduce and briefly discuss the / (i?, T^j gravitational 
models, where is the trace of the stress-energy of the scalar field. Some cosmological applications of this model 
are also presented. 

Since in the present model the covariant divergence of the stress-energy tensor is non-zero, the motion of massive 
test particles is non-geodesic, and an extra acceleration, due to the coupling between matter and geometry, is always 
present. The equations of motion of test particles are obtained from a variational principle. The same variational 
principle can be used to investigate the Newtonian limit of the model, and the expression of the extra-acceleration is 
also obtained. We use the precession of the perihelion of the planet Mercury to obtain a general constraint on the 
magnitude of the extra-acceleration. 

The present paper is organized as follows. The field equations of / (i?, T) gravity are derived in Section [TTJ Some 
particular cases of the model are considered in Section Mil The case of the scalar fields is discussed in Section IIV1 
and a Brans-Dicke type formulation of the model is obtained. The equations of motion of massive test particles are 
derived in Section[V] where the Newtonian limit of the model is also obtained and analyzed. We discuss and conclude 
our results in Section [VTl In the present paper we use the natural system of units with G = c = 1, so that the Einstein 
gravitational constant is defined as k 2 = 8tt. 



II. GRAVITATIONAL FIELD EQUATIONS OF / (R, T) GRAVITY 



We assume that the action for the modified theories of gravity takes the following form 



(1) 



where / (R, T) is an arbitrary function of the Ricci scalar, R, and of the trace T of the stress-energy tensor of the 
matter, T^ v . L m is the matter Lagrangian density, and we define the stress-energy tensor of matter as f26j 



2 HV-9 L m) 
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and its trace by T — g^T^, respectively. By assuming that the Lagrangian density L m of matter depends only on 
the metric tensor components <? M „, and not on its derivatives, we obtain 



9L m 



(3) 



By varying the action S of the gravitational field with respect to the metric tensor components 9 ^ v provides the 
following relationship 



ss 



1 

167T 



f R (R,T)SR + f T (R,T) -^5g^ -±g^f (R,T) 5<T + ltor- M v ! ' Lm) 



f=g Sg^ 



-gd 4 x, (4) 



where we have denoted f R (R, T) = df (R, T) /dR and f T (R, T) = df (R, T) /dT, respectively. For the variation of 
the Ricci scalar, we obtain 



5R = S (g^R^) = R^Sg^ + (V X ST^ - VJT 



(5) 



where Va is the covariant derivative with respect to the symmetric connection T associated to the metric g. The 
variation of the Christoffel symbols yields 



i 



g a (V^dgvu + V u 5g afi - V a Sg^ v ) , 



and the variation of the Ricci scalar provides the expression 

SR = R^SgT + Q^nSgT - V^V* . 
Therefore, for the variation of the action of the gravitational field we obtain 
1 



SS 



16tt 
+fr(R,T) 



f R (R, T) R^Sg^ + f R (R, T) g^Ubg^ - f R (R, T) V^/Jg^ 



S {g al3 T afj ) 1 1 8 



'—gd x . 



We define the variation of T with respect to the metric tensor as 

S {g aP T aP " 



TjXV ® fits ■ 



(6) 
(7) 



(8) 



(9) 



where 



= g 



t pSTa/3 



(10) 



After partially integrating the second and third terms in Eq. ©, we obtain the field equations of the / (R, T) gravity 
model as 

f R {R, T) R„ v - -/ (R, T) 9lu , + ( 9flv D - V^ u ) f R (R, T) = 8ttT^ - f T (R, T) - f T (R, T) . (11) 

Note that when f(R,T) = f(R), from Eqs. (jlll) we obtain the field equations of f(R) gravity. 

By contracting Eq. (JTTJ gives the following relation between the Ricci scalar R and the trace T of the stress-energy 
tensor, 



f R (R, T)R + 3Df R (R, T) - 2/ (R, T) = 8ttT - f T (R, T)T - f T (R, T) 6 : 



(12) 



where we have denoted 9 = 0A 

By eliminating the term Df R (R, T) between Eqs. (TTTj) and (fl^j) . the gravitational field equations can be written in 
the form 



1 



1 



f R (R, T) R pv - -Rg^ + -f {R, T) 9fiV = Sir - -Tg^ - f T (R, T) - -Tg^ 



1. 



1. 



-S T (r, t) [ e M „ - -e 9fiV ) + v M v„/ fl (r, t) 



(13) 
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Taking into account the covariant divergence of Eq. (TTT|) . with the use of the following mathematical identity [13] 

0, (14) 



f R (R, T) - If (R, T) giu> + (g^D - V M V, y ) f R (R, T) 



where / (i?, T) is an arbitrary function of the Ricci scalar R and of the trace of the stress-energy tensor T, we obtain 
for the divergence of the stress-energy tensor the equation 



h(R,T) 
8n-f T (R,T) 



V ^ = J f _)J^ [(Tp, + 6 M „) V In fx (R, T) + V6 F ] . (15) 



Next we consider the calculation of the tensor M „ , once the matter Lagrangian is known. From Eq. we obtain 
first 

STgfi _ Sg afi dL m d 2 L m 

6gi»> ~ 6gn» m + 9afi dg^ dg^dg°>P 

Sgafi T 1 1 d 2 L m 



From the condition g a ag' 7 ^ = 8^, we have 

6g a f 



9**9^512 , (17) 



SgV 

where 6^2 = Sg 17 " 1 /Sg^" is the generalized Kronecker symbol. Therefore for Q^ v we find 



e„„ = -2T, U + 9lxv L m - ^qP^ ■ (18) 

In the case of the electromagnetic field the matter Lagrangian is given by 

L m = --^F a0 F ja g^g^, (19) 

where F a p is the electromagnetic field tensor. In this case we obtain M „ = —T^ v . In the case of a massless scalar 
field (j) with Lagrangian L m = g a P\7 a <j>V p(j>, we obtain = —T^ v + (l/2)Tg^. The problem of the perfect fluids, 
described by an energy density p, pressure p and four-velocity is more complicated, since there is no unique 
definition of the matter Lagrangian. However, in the present study we assume that the stress-energy tensor of the 
matter is given by 

T^u = (p + p) u^u v - pg^ , (20) 

and the matter Lagrangian can be taken as L m = —p. The four-velocity satisfies the conditions u^u^ = 1 and 
u^VvU^ — 0, respectively. Then, with the use of Eq. (fl~8)). we obtain for the variation of the stress-energy of a perfect 
fluid the expression 

On» = -2T^ - pg^ . (21) 

III. PARTICULAR CASES OF GRAVITATIONAL FIELD EQUATIONS IN THE f(R,T) MODEL 

In the present Section we consider some particular classes of f(R, T) modified gravity models, obtained by explicitly 
specifying the functional form of /. Generally, the field equations also depend, through the tensor O^, on the physical 
nature of the matter field. Hence in the case of f(R,T) gravity, depending on the nature of the matter source, for 
each choice of / we can obtain several theoretical models, corresponding to different matter models. 

A. f(R,T) = R + 2f(T) 

As a first case of a f(R,T) modified gravity model we assume that the function f(R,T) is given by f (R,T) 
II + 2/(T), where f(T) is an arbitrary function of the trace of the stress-energy tensor of matter. The gravitational 
field equations immediately follow from Eq. ([lip, and are given by 

- \R9„u = 8ttT^ - 2f (T) T M „ - 2f\T)Q llu + f(T) 9fw , (22) 



5 



where the prime denotes a derivative with respect to the argument. 

If the matter source is a perfect fluid, M „ = —2T llv — pg^ vi then the field equations become 

- \R9»v = SttT^ + 2f (T) + [2pf'(T) + f(T)] g^ v . (23) 

In the case of dust with p — the gravitational field equations are given by 

R/J.V - -Rg^u = 8ttT^ + 2/'(T)T MV + / (T) 5ai! , . (24) 

These field equations were proposed in [25j to solve the cosmological constant problem. The simplest cosmological 
model can be obtained by assuming a dust universe (p — 0, T = p), and by choosing the function f(T) so that 
f(T) = AT, where A is a constant. By assuming that the metric of the universe is given by the flat Robertson- Walker 
metric, 

ds 2 = dt 2 - a 2 (t) (dx 2 + dy 2 + dz 2 ) , (25) 
the gravitational field equations are given by 

a 2 

3— = (8tt + 3A)p, (26) 

+ K = Ap, (27) 
a a z 

respectively. Thus this f(R,T) gravity model is equivalent to a cosmological model with an effective cosmological 
constant A e g oc H 2 , where H = d/a is the Hubble function 25j. It is also interesting to note that generally for 
this choice of f(R,T) the gravitational coupling becomes an effective and time dependent coupling, of the form 
G ff = G ± 2/'(T). Thus the term 2/(T) in the gravitational action modifies the gravitational interaction between 
matter and curvature, replacing G by a running gravitational coupling parameter. 
The field equations reduce to a single equation for H, 

8?r + 2A , 

2H + 3 ^T3X H =°> < M > 



with the general solution given by 



The scale factor evolves according to a(t) — t a , with a — 2 (87r + 3A) /3 (8ir + 2A). 



B. f(R,T) = f 1 (R) + f 2 (T) 

As a second example we consider the case in which the function / is given by / (R,T) = fi(R) + fiiT), where 
fi(R) and fe(T) are arbitrary functions of R and T, respectively. In this case for an arbitrary matter source the 
gravitational field equations are given by 

f[(R)R^ - \fi{R)9„ v + (g^a - V M V„)/i(J2) = B*T^ - ^(T)T^ - f 2 {T)Q^ + \hT)g^ . (30) 
Assuming that the matter content consists of a perfect fluid, the gravitational field equations become 



f[(R)R^ - \h{R)g^ + (g^n - v„v„)/{(fl) = 8nT^ + ftfryr^ 



f2(T)P+\f2(T) 



9^ 



(31) 



In the case of dust with p = 0, the gravitational field equations reduce to 



f[{R)Rpv - \h{R)g»» + (g^a - v„v u )f[(R) = SttT^ + &(t)t^ + l -h{T)g^ . 



(32) 
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In the case fi(T) = 0, we re-obtain the field equations of standard f(R) gravity. Eq. (|3Tj) can be reformulated as an 
effective Einstein field equations of the form 



G 



fits 



R 



1 



fJ,U 



where we have denoted 



m) 



8 71 



and 



(33) 



(34) 



(35) 



The gravitational coupling is again given by an effective, matter (and time) dependent coupling, which is proportional 
to the derivative of the function f% with respect to T. The gravitational field equations can be recast in such a form 
that the higher order corrections, coming both from the geometry, and from the matter-geometry coupling, provide 
a stress-energy tensor of geometrical and matter origin, describing an "effective" source term on the right hand side 
of the standard Einstein field equations. In the f(R,T) scenario, the cosmic acceleration may result not only from 
a geometrical contribution to the total cosmic energy density, but it is also dependent on the matter content of the 
universe, which provides new corrections to the Hilbert-Einstein Lagrangian via the matter-geometry coupling. 
The (t, t) component of Eq. (|5"!?|) has the following form: 



3H 2 = 



f'i(R) 



l 



8tt 



1 



fl(R) 



-~ (fi(R) -6(6+ 2H 2 ) f[{R)) + 2f' 2 {T) -)( /'/ • Mill) /['I /?.) 



(36) 



(37) 
(38) 



Here R = 6 yH + 2H 2 j . For simplicity, correspondsto the matter with a constant EoS parameter w. If we now 
define the e-folding TV by a = aoe N , p and T are given by 

p = p e- 3(1+tu)JV , T = -(1 - 3 W )poe- 3(1+tu)JV . 

We now consider an arbitrary development of the expansion in the Universe given by 

H = h(N) , 

where h(N) is an arbitrary function of N. Then Eq. (f36f can be written as 
fi(T) = F 2 (N) 

- 1 + poe 3 3(1+ ^ {- Y Poe~^ N + \h [6 (h(N)h'(N) + 2h(N) 2 )] 
- [h{N)h'{N) + h(N) 2 ] f[ [6 (h(N)h'(N) + 2h(N) 2 )] 

+3 [h(N) 2 h"(N) + h(N)h'(N) 2 + 4h(N) 2 h'(N)] f{ [6 (h(N)h'(N) + 2h(N) 2 )] } . 
Eq. (|39)l dictates that for an arbitrary fi(R), and for the following specific choice 



(39) 



&(T)=F 2 - 



In 



(1— 3w)p 



3(1 + w) 



(40) 



an arbitrary development of the expansion in the Universe given by (|38p can be realized. Hence, for viable f(R) 
gravitational models, using the above reconstruction method, the possibility arises to modify the universe evolution 
by adding the corresponding function depending on the trace of the stress-energy tensor. 



C. f(R,T) = MR) + f 2 (R)MT) 



As a third case of generalized f(R,T) gravity models, we consider that the action is given by / (R,T) = 
fi(R) + f2{R)}z{T), where /j, i = 1, 2, 3 are arbitrary functions of the argument. For an arbitrary matter source the 
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gravitational field equations are given by 

[fi(R) + f' 2 {R)h{T)\ - \h{R)9v» + {9^U - V M V„) [fi(R) + f 2 {R)h{T)\ 

= SttV - h{R)&{T)Tp, - f 2 (R)&(T)Q^ + \h{R)h{T) 9llv . (41) 
In the case of a perfect fluid we find the field equations 

[f[{R) + f^(R)MT)} R^ - \h{R)g^ + (g^D - V M V„) [f[(R) + f' 2 {R)h{T)] 



= SttT^ + f 2 (R)f^T)T^ + f 2 (R) 

For the case of dust matter we obtain 



&{T)P + \h(T) 



9n» ■ (42) 



{f[(R) + &(R)fa(X)] - \h{R)g^ + (g^D - V M V„) [f[(R) + f 2 (R)f 3 (T)} 

= SttT^ + f 2 (R)f^(T)T^ + \f 2 {R)h{T) 9llu . (43) 

In this class of models both, the effective cosmological constant A e g and the running gravitational coupling G e g are 
functions of both matter and geometry. 

IV. / (R, T 4 ") GRAVITY 

Scalar fields are supposed to play a fundamental role in physics and cosmology [28j |. In particular, cosmological 
inflation, the late-time acceleration of the universe, or dark matter and its properties can be explained in the framework 
of specific scalar field models. However, obtaining more general gravitational models with scalar fields as a source may 
give a better insight in the general properties of the gravitational field, and could also provide some possibilities for 
observationally testing the generalizations of gravity models. In the present Section, we consider the f(R,T) gravity 
model in the case of self-interacting scalar fields. 

A. The action of the / (R,^) gravity 
We start with the following action for matter, 

^matter (#/«/, Ipi) = / (fxy/^L (g pa , Ipi) . (44) 



In Eq. (j44|) the ipi's, i = 1,2, ... represent the matter fields. By using the matter action (|44[) . we now introduce the 
action for the gravitational field with matter sources as 

S = 2^ / d A x^g~F{R, dp) + 5 matt or (e+g^, , (45) 

where F(R, </>) is an algebraic function of R and of the scalar field dp. Then by the variation of the action with respect 
to dp, we obtain first 

^ F<l> (R,4>) + ln = 0, (46) 

where we denoted F^,(R, dp) = dF(, dp)/ddp and 

T EE e~V"%^ T 0/^ = T ^\g^^g^ > (4 7 ) 

respectively. The stress-energy tensor of matter is defined, as usual, by 
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By the assumption that F(R, <fi) is an algebraic function of R and <f>, Eq. (|46|) can be algebraically solved with respect 
to <j). Thus we can obtain <j> as a function of R and X^, i.e., (j) — <j) {R, T^). Then by substituting the expression of <p 
into the action (|45|) . we obtain an example of F(R,T^) gravity, with the following action 



S= 2^ I d x ^9 F ( R > T <l> 



(49) 



where we have denoted 

F (-R> T$) = F[R,<j> (R, Tcj,)} . (50) 
With the use of the conformal Weyl transformation g^ u -+ e _ ^<7 Mt/ , the action (|4"5)) or (14TH) is transformed as 



d 4 Xy/^g~F 



1 

2^ 



e*T 



^matter 
+ ^matter (<7/ii/, V'i) i 



(51) 



with T = g^T^y. In the action S'matter (<7/u>, i n Eq. (|51|) . the matter fields have only a minimal coupling with 
gravity, and they do not couple with <f>. Then the frame in the action (1511) might be regarded as a physical frame. 



B. Example of / (7?, T^) scalar field gravity, and reconstruction 



As an example of / (R, T^) gravity of the form R+f (T^) , we consider the case of a scalar field with a self-interaction 
potential V (</>)■ The action is given by 



S*= d A xJ=g~ 



--wicfyd^d^-V^) 



(52) 



where we have included u(4>) for later convenience. For the scalar field model described by Eq. (|52"1) . the trace of the 
stress-energy tensor is given by 



T* = -u{<j>)d^<l> - 4V(<f>) . 
Consequently we may define the / (i?, T^J — R + f (T^J gravity model in the following form: 



S 



d x\f— t 



(53) 



(54) 



For the model (|54|) . in a flat Friedman- Robertson- Walker geometry, the Friedman equations have the following form: 



\h 2 = i^(0)0 2 + y(0)-/L(0)0 2 -4y(0) 



2/' 



u{(j))4> 2 - W((j>) w(0)0 2 



-4 ( 3H Z + 2H 



1 

2 W 



V(<P) + f uo^ 2 -W{cj>) 



(55) 
(56) 



where H = a/a. 

In the following, we consider, for simplicity, the case V((f>) = 0. Then the action ([54]) has the following form: 



5 = | d i xJ—g^R + F[-u{4>)d^d^] 



and the Friedman equations (|55|) and (|56p take the form, 

-F 



2F' 



u;(0)0 2 u>(#/> 2 , 



1 



3# 2 + 2iJ 



F 



w(</# 5 



(57) 

(58) 
(59) 
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In the action (l57l) . F [— 0^(0)9^09^0] is denned by 



F = / [-a;(0)9 M 0^0] - ^(0)^0^*0 . (60) 



The action (|57|) gives a model of k-essence |29M31| . In [3^j , it has been shown that the Friedmann equations (|55|) and 
([59| do not admit the de Sitter solution, except in the trivial case where is a constant, and F(0) > 0. In [32|, the 
formalism of the general reconstruction has also been explicitly given. An explicit model of modified gravity in which 
a crossing of the phantom divide can be realized was reconstructed in (33j . 
As a simple example, we consider the model 

F [-oj(0)^0^0] = -F Q e- 2ln (^) d ^ d ^ , (61) 
where F$ and 0o are constants. The Friedman equations have a solution where the universe expands by a power law, 

H = Y> <t> = L (62) 

The constant ho can be obtained by solving the following algebraic equation 

3hl - 2h = k 2 0o^o • (63) 

V. THE EQUATION OF MOTION OF TEST PARTICLES AND THE NEWTONIAN LIMIT IN / (R, T) 

GRAVITY 

Since in the general f(R,T) type gravity models the stress-energy tensor of matter is not covariantly conserved, it 
follows that the test particles, moving in a gravitational field, do not follow geodesic lines. This situation is similar 
to the case of the / (R,L m ) models [24|, where the coupling between matter and geometry induces a supplementary 
acceleration acting on the particle. In the present Section, we derive the equations of motion of test particles in 
f(R,T) gravity models, and obtain the Newtonian limit of the theory. We also investigate the constraints on the 
magnitude of the extra-acceleration that can be obtained from the available observational data on the perihelion 
precession of the planet Mercury. 

A. The equations of motion of test particles 



In the case of a perfect fluid, with the stress-energy tensor given by Eq. (|20l) . the divergence of the stress-energy 
tensor is given by 

VT M „ = - / {T^f T (R,T)+g^[f T (R,T)p]} . (64) 

ott + It (R, L ) 

We also introduce the projection operator — g„\ — u^u\ for which we have h^xu 11 = and h fl \T ilv = —h^p, 
respectively. 

Explicitly, Eq. (|64p can be written in the form 

V, (p + p) u^u" + (p + P ) [v?V v u» + u^ v u v ] - g^V„p 

= - « 177^ {T^Vuh (R, T) + gTV v [f T (R, T)p]}. (65) 

OTT + JT (R, L ) 

By contracting Eq. (|65l) with we obtain 

After multiplying with g aX and by taking into account the identity 
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we obtain the equation of motion of a test fluid in / (R, T) gravity as 

rO'V = F , (68) 



where 



ds 2 

f " = 877 ( Z U« V ff (n TM ^ ~ U ^ UU) ■ (69) 
(p + p) [8tt + f T (R, T)\ 

The extra-force / M is perpendicular to the four-velocity, f^u^ = 0. When fx (R, T) = 0, we re-obtain the equation of 
motion of perfect fluids with pressure in standard general relativity, which follow from the conservation of the energy- 
momentum tensor, V = 34]. In the limit p 0, corresponding to a pressureless fluid (dust), in standard 
general relativity the motion of the test particles becomes geodesic. The same result holds true in the f(R, T) gravity 
model. Even if fx (R, T) ^ 0, the motion of the dust particles always follows the geodesic lines of the geometry. By 
assuming that the term %irV v p/ (p + p) [8ir + fx (R, T)] can be formally represented as V„ In \fQ, 

8tt , - Ug V ^ f m =V„lnVQ, (70) 
(p + p) [8tt + f T (R, T)\ 

the equation of motion Eq. (|68[) can be obtained from the variational principle 



SS P = S J L p ds = 6 ^Qy/g^Wda = , (71) 

where S p and L p = ^/Q^/g^u^u 1 ^ are the action and the Lagrangian density for the test particles, respectively. 
To prove this result we start with the Lagrange equations corresponding to the action (|T1 [) , 

d_ ( dLp\ _ dLp _ 
ds \du x J dx x 

Since 



dL p 
du x 



(72) 

= (73) 



and 



a straightforward calculation gives the equations of motion of the particle as 

+ + {u»u v - g^) V„ln y/Q = . (75) 

When y[Q — > 1 we re-obtain the standard general relativistic equation for geodesic motion. 

As an example of the application of the previous formalism we consider the case in which the pressure can be 
expressed as a function of the density by a linear barotropic equation of state of the form p = wp, where the constant 
w satisfies the condition w -C 1. Therefore p + p rj p and T = p — 3p w p, respectively. Moreover, for simplicity, 
we also assume that the function fx is a function of T p only. We can expand fx near a fixed value po of the 
density, so that fx (ft) = fx (fto) + (ft- Po) [dfx /dp) /dp] \ p=Po = 8ir [o + b (p - po)], where a = fx (po) /Stt and 
bo = [dfx I dp) /dp] |p =Po /87r, respectively. Eq. (f70| of the definition of y/Q becomes 

; - W t Vvln — / h r = V„h^, (76) 

1 + ao - o po 1 + ao + bo (p - po) 

giving 



(P) 



Cp 



-| w/(l+a -b p ) 

(77) 



1 + a + bo (p - po)_ 

where C is an arbitrary constant of integration. Eq. (|70|) is also valid for a fluid satisfying a linear barotropic equation 
of state of the form p = (7— 1) p, 7 = constant, and for a model with fx(R,T) = constant = fx- In this case 
y/Q = C , T/0 87r< - 7_1 - | / 7( - 87r+ ^ T -', where Cx is an arbitrary integration constant. Therefore Eq. (1701) is valid in both the 
non-relativistic and the extreme relativistic limits of the model. On the other hand we have to mention that the 
function ^/Q can always be obtained by formally integrating the left-hand side of Eq. (I70|) . However, generally this 
function cannot be expressed in an exact analytical form, and to find its functional form approximate methods have 
to be used. 
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B. The Newtonian limit 



The variational principle (fTTj) and the pressureless dust model, described by Eqs. ([76]) and (l77l) . can be used to 
study the Newtonian limit of the model. In the limit of the weak gravitational fields, 



ds w a/1 + 20 — &dt w (1 + - t/ 2 /2) di . 



(78) 



where is the Newtonian potential and v is the usual tridimensional velocity of the fluid. By using the relation 
x a — exp(alna;) ~ 1 + ahix, we can approximate \[Q (p) given by Eq. (1771) as 



Q(p)wl + 



(1 + oo - 6oPo) 



In 



Cp 



where we have denoted 



c/(p) 



1 + a + b (p - po) 
Cp 



1 + U(p) 



(1 + a - b p ) 



hi 



1 + a + b a (p - po) 



(79) 



(80) 



In the first order of approximation the equations of motion of the fluid can be derived from the variational principle 



i + u ( P ) + <t> - 



dt = 0. 



ID 



and are given by 

a = -V(j> - VJ7 (p) = aw + d p + d E , (82) 
where a is the total acceleration of the system, ajv = — V0 is the Newtonian gravitational acceleration and 



C 1„ 1„ 

— Vp = Vp, 

1 + ao - bopo p p 



(83) 



is the hydrodynamical acceleration. Eq. (|53")) also allows us to fix the value of the arbitrary integration constant C as 
C = 1 + ao — 6qPo- Finally, 



as (p,p) = 



Vp 



1 + a - b po l + ao + b (p- p ) ' 
is a supplementary acceleration induced due to the modification of the action of the gravitational field. 



(84) 



C. The precession of the perihelion of Mercury 

An estimation of the effect of the extra- force, generated by the coupling between matter and geometry, on the orbital 
parameters of the motion of the planets around the Sun can be obtained in a simple way by using the properties of 
the Runge-Lenz vector, defined as A = v xL — ae r , where v is the velocity relative to the Sun, with mass Mq, of a 
planet of mass m, r = re r is the two-body position vector, p = pv is the relative momentum, p = tuMq/ (to + Mq) 
is the reduced mass, L — f xp — pr 2 9k is the angular momentum, and a = GuiMq (35j | . For an elliptical orbit of 
eccentricity e, major semi-axis a, and period T, the equation of the orbit is given by (L 2 /pa) r" 1 = 1 + ecosO. The 
Runge-Lenz vector can be expressed as 



A 



a e r — rLeg . 



(85) 



and its derivative with respect to the polar angle 9 is given by 



dA 
~dJ 



~dV{r) 
dr 



a 



('o ■ 



(86) 
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where V(r) is the potential of the central force [Hj]. The potential term consists of the Post-Newtonian potential, 
Vpw(r) — ~a/r — 3a 2 /mr 2 , plus the contribution from the general coupling between matter and geometry. Thus we 
have 



dA 



ma E {r) 



(87) 



where we have also assumed that /i rs m. The change in direction Acf> of the perihelion with a change of 9 of 2ir is 



obtained as A4> = (1/ae) f*Lx dA/d6 



d9, and it is given by 



= 24tt ; 



i (i 



.3/2 



e 2 87r 3 me (o/T) Jo 



2tt ag 



L 2 (l + ecos#) 1 /ma 



(1 



cos 0cf0 . 



(88) 



where we have used the relation a/L — 2ir (a/T) /yl — e 2 . The first term of this equation corresponds to the standard 
general relativistic precession of the perihelion of the planets, while the second term gives the contribution to the 
perihelion precession due to the presence of the coupling between matter and geometry. 

As an example of the application of Eq. (|551) we consider the case for which the extra-force may be considered as 
a constant, cie ~ constant, an approximation that could be valid for small regions of spacetime. In the Newtonian 
limit the extra-acceleration generated by the coupling between matter and geometry can be expressed in a similar 
form [l!| . With the use of Eq. (1551) one finds for the perihelion precession the expression 



Ac/) 



6-irGM, 



G 



27ra 2 vT 



a(l-e 2 ) 



GMr. 



-CLE ■ 



(89) 



and 



where we have also used Kepler's third law, T 2 = 4tt 2 a 3 /GM Q . For the planet Mercury a = 57.91 x 10 11 cm, 
e = 0.205615, respectively, while M Q = 1.989 x 10 33 g. With these numerical values the first term in Eq. 
gives the standard general relativistic value for the precession angle, (A(j)) GR = 42.962 arcsec per century, while 
the observed value of the precession is (A0) obs = 43.11 ± 0.21 arcsec per century [36| . Therefore the difference 
(A(j)) E = (A<fi) obs — (A<fi) GR — 0.17 arcsec per century can be attributed to other physical effects. Hence the 
observational constraints requires that the value of the constant a E must satisfy the condition a E < 1.28 x 10~ 9 
cm/s 2 . 



VI. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have considered a generalized gravity model with an arbitrary coupling between matter 
(described by the trace of the stress-energy tensor) and geometry, with the Lagrangian given by an arbitrary function 
of T and of the Ricci scalar. We have derived the gravitational field equations corresponding to this model, and 
considered several particular cases that may be relevant in explaining some of the open problems of cosmology and 
astrophysics. The new matter and time dependent terms in the gravitational field equations play the role of an effective 
cosmological constant. We have also demonstrated the possibility of reconstruction of arbitrary FPW cosmologies by 
an appropriate choice of a function f(T). The equations of motion corresponding to this model show the presence 
of an extra-force acting on test particles, and the motion is generally non-geodesic. We have obtained, by using the 
perihelion precession of Mercury, an upper limit on the magnitude of the extra-acceleration in the Solar System. This 
value of aE, obtained from the solar system observations, is somewhat smaller than the value of the extra-acceleration 
a E Ri 10~ 8 cm/s 2 , necessary to explain the "dark matter" properties, as well as the Pioneer anomaly [19l. I37L l38j|. 
However, it does not rule out the possibility of the presence of some extra gravitational effects acting at both the solar 
system and galactic levels, since the assumption of a constant extra-force may not be correct on larger astronomical 
scales. 

Therefore the predictions of the f(R,T) gravity model could lead to some major differences, as compared to the 
predictions of standard general relativity, or other generalized gravity models, in several problems of current interest, 
such as cosmology, gravitational collapse or the generation of gravitational waves. The study of these phenomena 
may also provide some specific signatures and effects, which could distinguish and discriminate between the various 
gravitational models. In order to explore in more detail the connections between the f(R,T) gravity model and the 
cosmological evolution, some explicit physical models are necessary to be built. This will be done in forthcoming 
work. 
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